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INTRODUCTION 


There  is  an  extensive  literature,  spanning  more  than  four  decades,  des¬ 
cribing  investigations  of  spatially  localized  streams  or  plasma-magnetic 
field  entities  (plasmoids)  in  interaction  with  background  magnetic  field 
environments.  Those  investigations  include  theoretical  analyses,  com¬ 
puter  simulations,  and  experimental  research  efforts.  The  large  fraction 
of  these  investigations,  however,  have  concerned  plasma  interactions 
characterized  by  MA<1,  where  is  the  Alfven  mach  number  defined 
such  that: 

2  2 
*  (plasma  flow  velocity/Alfven  speed) 

=  (plasma  kinetic  energy  density) /(magnetic  field  energy 
density) 

It  is  well  established  that  in  such  regimes,  the  plasma  is  either 
strongly  decelerated  (and/or  trapped)  by  the  transverse  magnetic  field  or 
propagates  unimpeded  (only  after  undergoing  charge  polarization  effects) 
at  the  reduced  ExB  drift  velocity.  Very  little  work  has  been  done  in  the 
"super-Alfvenic"  regimes:  M.»l.  This  is  the  regime  of  critical 

A 

importance  in  the  exoatmospheric  propagation  of  very  energetic  neutral¬ 
ized  charged  particle  beams.  Recent  computer  simulations  in  that  regime 
(e.g. ,  those  of  the  theory  group  effort  at  Los  Alamos  National  Laboratory 
(LANL))  have  shed  some  light  on  propagation  dynamics  but  have  been  handi¬ 
capped  by  the  following  limitations: 
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a.  The  lack  of  any  clearly  justifiable  choice  criteria  for  selecting 
among  the  many  possibilities  for  collective  mode  instabilities  in 
that  growth  rates  depend  critically  upon  assumed  initial  plasma 
flow  structures. 

b.  The  computational  complexity  and  time  intensivity  of  full  three- 
dimensional  simulations  necessary  tc  treat  correctly  the  boundary 
effects  which  are  inescapable  in  actual  beam  propagation. 

c.  The  limited  ability  to  follow  the  plasmoid  for  long  (simulation) 
propagation  distances  (related  to  computer  memory  size/computa¬ 
tion  time  limitations)  and/or  the  inability  to  choose  rationally 
among  possible  computation  box  boundary  conditions. 

d.  An  absence  of  knowledge  of  the  significance  of  nonlinear  proc¬ 
esses  in  either  saturating  linearly  growing  instabilities  or 
aggravating  their  effect. 

Associated  theoretical  efforts  have  almost  always  been  restricted  to 
linearized  treatments  and/or  have  assumed  (for  tractability)  highly 
simplified  geometries. 


* 
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It  that  light,  a  carefully  organized  set  of  scaled  simulation  experiments 
can  thus  provide  important  new  knowledge  which  can  support  and  complement 
the  theoretical  and  computer  simulation  efforts.  This  is  the  signifi¬ 
cance  of  the  effort  supported  by  the  contract  (F49620-81-C-0021)  now 
ended.  The  experimental  facility  which  was  used  is  located  at  Morehouse 
College,  Atlanta,  Georgia  and  was  made  available  to  AMAF  Industries,  Inc. 
on  a  leased  basis. 


3  (and  4) 


STATEMENT  OF  WORK 


The  research  objectives  for  the  effort  as  proposed  were  defined  by  the 
following  tasks  (see  Appendix  A  for  the  Statement  of  Work  as  it  actually 
appeared  in  the  contract  document): 

1.  Establish  effective  channel  of  communications  with  existing  computer 
simulation  groups  (e.g.,  LASL) . 

2.  Review  published  literature  on  experimental  and  theoretical  research 
on  instabilities  in  the  plasma  flow  regime  M^>1  for  applicable  insights. 

3.  Reconfigure  existing  hypervelocity  plasma  flow  generator  facility  to 
operate  with  high  reproducibility  in  the  regime  MA»1  as  well  as 

This  will  involve  (as  simple  engineering  tasks): 

a.  Extension  of  the  vacuum  chamber  in  the  flow  direction  by  addition 
of  a  four-foot  pyrex  glass  pipe  section. 

b.  Redesign  and  reconstruction  of  the  magnetic  field  producing  coil 
system  (currently  a  Helmholtz  pair)  to  produce  an  approximately 
uniform  field  over  the  extended  flow  channel. 
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c.  Addition  of  a  ballast  chamber,  a  precision  leak  valve,  a  mechan¬ 
ical  vacuum  pump,  and  a  vacuum  gauge  to  supply  controllable  quan¬ 
tities  of  selected  gases  to  the  breech  of  the  gun  (operating  in 
the  "puff"  mode). 


4.  Construct  and  calibrate  fast  response  magnetic,  electric  and  current 
probes  for  flow  diagnostics  (designed  by  well-known  techniques,  to  have 
response  time  <,0.1/tsec). 


5.  Procure  and  bri.ig  to  operational  status  a  fast  optical  system  (image 
convert-*-  or  Kerj  optics)  for  macroscopic  structure  observations. 


6.  Begin  program  of  experimental  parmetric  surveys  (with  full  diagnostic 
monitoring) : 


a.  Plasma  beam  velocity. 


b.  Plasma  beam  density. 


c.  Plasma  beam  particle  mass  (variations  in  composition), 


d.  Propagation  chamber  ambient  pressure. 


Background  magnetic  field  intensity. 


f.  Background  magnetic  field  structure. 
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7.  Analyze  results  of  experimental  surveys  as  to: 


a.  Dependence  of  propagation  stability  on  parameter  values. 


b.  Characteristic  of  instability  structure  (growth  rates,  spatial 
scales,  onset  conditions,  saturation  levels). 


8.  Summarize  results  in  form  for  inclusion  in  computer  simulations 
(especially  in  computer  simulation  areas  a.  and  d.  described  previously). 


9.  Investigate  experimentally  (guided  by  appropriately  modified  LASL 
computer  code  predictions  as  they  become  available)  possible  means  for 
stabilizing  the  propagating  beam  with  an  aim  of  achieving  the  highest 
energy  ir  the  beam  consistent  with  stability. 


7  (and  8) 


STATUS  OF  THE  RESEARCH  EFFORT 


THE  FACILITY  CONVERSION: 


The  facility  utilized  was  located  in  the  Department  of  Physics  at  More¬ 
house  College  in  Altanta,  Georgia.  It  had  previously  been  used  in  a 
Department  of  Energy  sponsored  research  effort  in  pulsed,  inductive 
direct  energy  conversion  based  on  radial  expansion  flows.  The  major 
components  of  the  system  are: 

1.  The  coaxial  plasma  gun  and  associated  capacitor  bank. 

2.  The  magnetic  field  coil  system  and  associated  capacitor  bank. 

3.  The  vacuum  chamber  and  associated  pumps. 

4.  Plasma  and  field  diagnostics. 

5.  The  data  acquisition/control  room. 

A  large  area  view  of  the  above-table-top-visible  psrtion  of  the  facility 
is  provided  in  Figures  la  and  lb. 
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1.  The  Coaxial  Gun  and  Capacitor  Bank 

The  coaxial  plasma  gun  is  a  standard  Marshall  gun  configuration  capable 
of  being  either  metallic  foil  or  gas  "puff"  breech  loaded.  In  the  gas 
"puff"  mode  it  is  loaded  by  a  thyratron  controlled,  solenoid  actuated  gas 
valve.  The  gun  is  constructed  of  copper  outer  and  brass  rod  inner  elec¬ 
trodes  with  an  outer  radius/inner  radius  of  4  and  a  length  of  10  cm.  The 
plasma  composition  options  include  He,  N,  A,  or  various  metals  (from 
foils).  Details  of  the  construction  of  the  gun  is  provided  in  Figures  2a 
and  2b.  The  capacitor  bank  for  the  gun  is  comprised  of  8  capacitors 
rated  at  15  uf  and  20  KV  in  parallel  yielding  a  maximum  storage  of  24 
kilojoules.  The  energy  in  the  bank  in  switched  by  4  type  GL-7703  igni- 
trons  in  parallel  (Figure  3).  The.  capacitor  bank  switched  into  the  gun 
has  a  "ring"  frequency  of  approximately  10  kHz.  A  typical  storage  oscil¬ 
loscope  trace  of  the  gun  voltage  is  shown  in  Figure  4.  Nominal  plasma 

22  -3 

densities  and  flow  velocities  are  respectively  10  m  and  10  km/sec. 
Figure  5  is  an  open  shutter  photograph  of  the  beam. 


2.  The  Magnetic  Field  Coil  System  and  Capacitor  Bank 
The  transverse  magnetic  field  is  produced  by  a  rectangular  coil  array 
extending  most  of  the  length  of  the  vacuum  chamber.  The  coils  are 
visible  in  Figures  la  and  lb.  The  coils  consist  of  4 'turns  each  and  are 
145  cm  long,  30  cm  wide  and  are  separated  by  25  cm.  The  array  has  an 
effective  inductance  of  approximately  100  uh.  The  array  is  powered  by  a 
capacitor  bank  comprise  of  4  capacitors  in  parallel  each  rated  at  60  uf 
and  20  KV.  The  bank  is  switched  by  4  type  GL-7703  ignitrons  in  parallel. 
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The  calculated  field  profiles  (normalized  to  the  field  at  the  array 
center)  are  plotted  in  Figures  6a  and  6b  for  the  orientation:  x,  width; 
y,  length;  z,  separation.  Figure  7  shows  the  array  location  on  the 
vacuum  chamber.  The  capacitor  bank  "dumps"  into  the  array  with  a  pulse 
(critically  damped)  that  peaks  in  200  usee  and  zeroes  in  500  usee. 

Figure  8  shows  a  typical  field  signal  as  picked  up  by  a  magnetic  field 
probe  at  the  array  center.  With  3  KV  on  the  bank,  the  field  at  the 
center  is  measured  to  be  approximately  0.05  weber/m  (500  gauss).  The 
fieid  is  proportional  to  the  bank  charge  voltage  over  a  large  range  of 
initial  charge.  An  estimate  of  the  fields  required  to  operate  in  the 
Ma>?1  regime  is  calculated  in  Appendix  D. 

3.  The  Vacuum  Chamber  and  Pumps 

The  vacuum  chamber  is  constructed  out  of  6"  diameter  Corning  pyrex  sani¬ 
tary  drain  pipe  sections  with  rubber  gaskets.  Figure  9  shows  the  overall 
dimensions  of  the  chamber.  The  gun  is  mounted  on  the  axis  of  the  "T"  end. 
The  chamber  is  pumped  by  a  large,  two-stage  mechanical  pump  (Welsh  1397) 
and  a  small  oil  diffusion  pamp  with  backing  pump  to  the  1-10  uHg  range. 

4.  Plasma  and  Field  Diagnostics 

The  background  and  short  time-scale  magnetic  fields  are  measured  by 
magnetic  probes  connected  to  passive  RC  integrators.  Figure  1C  shows 
such  a  probe  inserted  through  ?n  access  port  into  a  crossed-T  section.  A 
typical  calibration  for  such  a  probe  working  into  a  100  usee  RC  inte- 
grator  is  0.5  mV /gauss  (5  V/weber/m  ). 


The  plarma  velocity  is  measured  by  "time-of-f light"  techniques  using  a 
bifurcated  light  pipe  and  a  very  fast  response  photodiode  detector 
(Tropel).  Figure  11  shows  the  setup  and  Figure  12  shows  the  location  of 
the  velocity  probe  as  well  as  the  field  probe  on  the  chamber. 


The  density  gradient  structure  of  the  beam  is  measured  by  a  laser 
schliren  system  (Kiefer  &  Lutz,  1965)  with  a  2  m  travel  path  and  using  a 
5  mwatt  helium-neon  laser  as  a  source.  *'he  geometry  of  the  setup  is 
shown  in  Figure  13  and  portions  of  the  actual  system  are  shown  in  Figures 
14a,  14b,  and  14c.  In  this  system  the  voltage  from  the  photodiode 
detector  is  proportional  to  the  plasma  (local)  density  gradient. 


A  fast  optical  Kerr  cel’  shutter  system  has  been  designed  and  constructed 
but  is  not  yet  operative.  This  system  will  be  used  to  obtain  very  short 
time  interval  photographs  (100  nsec)  which  will  catch  the  beam  at  early 
times  in  its  propagation  and  allow  for  analysis  of  the  mOv.e  structure  of 
any  dynamic  instabilities. 


5.  The  Data  Acquisition/ Control  Room 

The  control  room  is  completely  enclosed  by  a  double  screen  shielded  room. 
Figure  15  shows  the  back  of  that  enclosure  and  the  BNC  data/signal  ports. 
The  control  room  instrumentation  includes  4  Tektronix  storage  oscillo¬ 
scopes  (Figure  16a)  for  data  monitoring.  Precisely  controllable  trigger 
signals  are  generated  by  the  combination  of  a  Tektronix  161/162  signal 
generator  and  an  Abtronics  Model  100  4  channel  time  delay  generator 
(Figure  16b). 
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Figure  6a 


Figure  6b 
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Figure  7 
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Knife  Edge 


Phocodiode 


Figure  15 
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STATUS  OF  THE  EXPERIMENTAL  EFFORT: 


Tasks  1-4  of  Che  Work  Statement  have  been  completed  as  projected.  Task  5 
(re  the  fast  optical  system)  is  still  in  a  design  implementation  and 
testing  phase.  Task  6,  the  experimental  parametric  surveys,  has  just 
begun.  Preliminary  results  appear  to  show  propagation  with  little 
velocity  attenuation,  beam  disruption,  or  magnetic  field  perturbation. 
This  mode  of  propagation  would  be  consistent  with  the  analysis  of  Baker 
and  Hammel  (D.  A.  Baker  and  J.  E.  Harnmel,  "Experimental  Studies  of  the 
Penetration  of  a  Plasma  Stream  in  a  Transverse  Magnetic  Field,"  Physics 
of  Fluids  8_,  713  (1965))  where  E  x  B  motion  plays  a  dominant  role  in 
propagation  dynamics.  It,  however,  is  important  to  emphasize  here  that 
these  results  are  only  preliminary.  A  careful  measurement  of  fast  time 
scale  magnetic  processes  during  the  field-flow  intermixing  region  of  the 
dynamics  of  propagation  has  yet  to  be  done.  Indeed,  the  first  year  of 
this  effort  has  been  necessarily  expended  in  bringing  the  facility  to  the 
operational  status  where  now  such  measurements  can  be  done. 


THE  THEORETICAL  EFFORT: 


Although  a  great  deal  of  theoretical  analysis  has  been  published  on 
across-field  propagation  of  plasma  beams,  the  large  majority  of  it  has 
begun  with  assumptions  (e.g.,  assumptions  of  polarized  E  x  B  steady 
states)  which  are  precisely  those  that  are  of  issue  here.  The  theo¬ 
retical  work  that  we  have  found  the  most  relevant:  W.  R.  Shanahan, 


"Energetic  Neutralized  Plasma  Beam  Propagation  Across  A  Magnetic  Field," 
Los  Alamos  Scientific  Laboratory,  LA-UR  81-2980,  itself  is  significantly 
limited  by  simplifying  assumptions.  Moreover,  the  assumptions  made  limit 
the  degree  to  which  the  analytical  solutions  derived  therein  could  be 
valid.  Indeed,  we  have  not,  even  after  careful  analysis  of  the  model  for 
the  propagation  there  presented,  been  able  to  come  to  Che  same  mathemati¬ 
cal  results.  We  have  initiated  a  dialogue  with  Dr.  Shanahan  to  resolve 
these  issues  and  to  advance  the  quality  of  the  theory.  We  propose,  and 
have  initiated,  a  more  adequate  theoretical  analysis  which  takes  the 
following  into  account: 

1.  The  early  time  competition  of  diamagnetic  and  polarization 
(dielectric)  response  processes  without  making  gross  quasi¬ 
steady  assumptions  (a  simple  transverse  dielectric  constant, 
for  example) . 

2.  The  possibility  of  large  amplitude,  nonlinear  evolution  of  the 
field. 

3.  The  critical  importance  of  beam  boundaries  on  beam  propagation 
dynamics. 

Such  theoretical  analysis  will  be  critical  to  any  experimental  effort  for 
interpretation  of  the  data  and  for  projecting  corrective  mechanisms 
should  catastrophic  instabilities  be  uncovered. 
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LIST  OF  JOURNAL  PUBLICATIONS 

There  have  been  no  journal  publications  resulting  from  the  research  as 
yet.  It  is  anticipated  that  results  of  interest  to  the  plasma  physics 
community  from  both  the  experimental  and  theoretical  efforts  will  be 
reported  at  the  American  Physical  Society  Plasma  Divisional  Meeting  in 
November  of  1982. 

PROFESSIONAL  PERSONNEL 

Principal  Investigator  -  Dr.  Carl  Spight 

Responsible  for  overall  research  direction,  technical  validity,  and 
project  management. 

Member  of  Technical  Staff  -  Dr.  Ronald  Graves 

Initially  responsible  for  computer  aided  design  of  magnetic  field  coils 
and  theoretical  analysis. 

Member  of  Technical  Staff  -  Dr.  Carlos  Handy 

Responsible  for  maintaining  theoretical  interaction  with  LANL  and  for 
pursuit  of  analysis  of  bounded  plasma  propagation. 

Member  of  Technical  Staff  -  Alfredo  Monge 

Responsible  for  execution  of  experiment  in  the  plasma  generator  facility 
at  Morehouse  College,  Atlanta,  Georgia. 
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INTERACTIONS  AND  PRESENTATIONS 

An  effective  channel  for  interaction  in  the  theoretical  effort  has  been 
established  at  Los  Alamos  National  Laboratory  (LANL)  through  Dr.  William 
R.  Shanahan  of  the  Applied  Theoretical  Physics  Division.  That  inter¬ 
action  has  allowed  us  to  clarify  the  strengths  and  limitations  of  his 
theoretical  and  numerical  simulation  analyses,  to  identify  assumptions  in 
that  work  of  a  problemmatic  nature,  and  to  define  a  reasonable  direction 
for  further  theoretical  analysis.  Analysis  of  the  stability  of  a  bounded 
plasma  beam  with  the  inclusion  of  polarization  effects  has  begun  based  on 
the  insights  gained  through  that  interaction. 

PROJECTIONS  FOR  FUTURE  EFFORTS 

Although  the  conversion  of  the  hypervelocity  plasma  generator  facility 
has  been  completed  and  a  full  set  of  diagnostics  implemented,  the  year  of 
effort  has  not,  however,  allowed  time  enough  for  meaningful  data  to  be 
gathered  on  the  stability  properties  of  the  beam  in  propagation  across 
the  applied  field.  Thus,  a  survey  of  beam  dynamics  as  a  function  of  gun 
and  field  parameters  remains  to  be  done.  In  addition,  the  absence  of  an 
adequate  theory  or  numerical  simulation  of  the  fully  bounded  beam  propa¬ 
gation  regime  points  to  the  importance  of  bringing  to  closure  the  theo¬ 
retical  effort  now  underway.  This  will  ensure  that  an  analytic  context 
will  be  available  against  which  to  make  sense  of  the  experimental  data. 
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PART  I  -  THE  SCHEDULE 
SECTION  B  -  SUPPLIES/SERVICES  AND  PRICES 


0001  RESEARCH 

The  contractor  shall  furnish  the  level  of  effort  specified  in  Section  F, 
together  with  all  related  services,  facilities,  supplies  and  materials  needed  to 
conduct  the  research  described  below.  The  research  shall  be  conducted  during  the 
period  specified  in  Section  F. 

0001AA 

Cl)  Establish  effective  channel  of  communications  with  existing  computer 
simulation  groups,  such  as  those  at  LANSL. 

(2)  Review  published  literature  on  experimental  and  theoretical  research 
on  instabilities  in  the  plasma  flow  regime  1  for  applicable  insights. 

(3)  Reconfigure  Morehouse  College  hypervelocity  plasma  flow  generator 
facility  to  operate  with  high  reproducibility  in  the  regime  majp^i  as  well 
as  maZ1. 

This  will  involve  (as  simple  engineering  tasks): 

(a)  Extension  of  the  vacuum  chamber  in  the  flow  direction  by  addition 
of  a  four  foot  pyrex  glass  pipe  section. 

(b)  Redesign  and  reconstruction  of  the  magnetic  field  producing  coil 
system  (currently  a  Helmholtz  pair)  to  produce  an  approximately 
uniform  field  over  the  extended  flow  channel. 


(c)  Addition  of  a  ballast  chamber,  a  precision  leak  valve,  a 

mechanical  vacuum  pump,  and  a  vacuum  gauge  to  supply  controllable 
quantities  of  selected  gases  to  the  breech  of  the  gun  (operating 
in  the  gas  "puff"  mode). 

(4)  Construct  and  calibrate  fast  response  magnetic,  electric  and  current 
probes  for  flow  diagnostics  (designed  by  well  known  techniques,  to  have 
response  times  0.1  sec.). 

(5)  Perform  experimental  parmetric  surveys '(with  full  monitoring)  of: 

(a)  Plasma  beam  velocity. 

(b)  Plasma  beam  density. 

(c)  Plasma  beam  particle  mass  (variations  in  composition). 
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d  Propagation  chamber  ambient  pressure, 

e  Background  magnetic  field  intensity, 

f  Background  magnetic  field  structure. 

6  Analyze  results  of  experimental  surveys  as  to: 

a  Dependence  of  propogation  stability  on  parameter  values. 

b  Characteristic  of  instability  structure  (growth  rates,  spatial 
scales,  onset  conditions,  saturation  levels). 

7  Summarize  results  in  form  for  inclusion  in  computer  simulations 
(especially  in  computer  simulation  areas  a.  and  a.  described  previously). 

8  Investigate  experimentally  (guided  by  appropriately  modified  computer 
code  predictions  as  they  become  available)  possible  means  for  stabilizing  the 
propagating  beam  with  an  aim  of  achieving  the  highest  energy  in  the  beam 
consistent  with  stability. 


Page  3  of  12  Pages 


AMAF,  INDUSTRIES,  INC. 
F49620-81-C-0021 


>  LIS 


CL  FIELD  13:47ED7  05/10/82 

100*  current  loop  magnetic  field  calculation 

I  110* 

120  integer  n»nx»ny,nz 

130  REAL  SVECT (5,7) » XHINf XHAX » YHINt YHftX rZMINf ZHftX » DX r DY *DZt h (5) .NORM 
140  OPEN (FILE3 ’ PTFIL’ ?  UN I T= i , STATUS =' UNKNOWN ’ > ACCESS-' LINENO’) 

150  CALL  SCRSUP (SVECT » H,H> 

160  CALL  GRIBUHIN, YMIN, XMAX, YMAXrZMIN, ZMAX, NX,NY,NZ»DX,DY,DZ) 

170  NORM=l ,0 

•  1S0  CALL  DFCALC (SVECT, M, 0,0, 0,0, 0,0, 1 . » 1 . » 1 . ,NX,NY,NZ,N»NQRM) 

190  CALL  DFCALC (SVECT  , M  >  XMIN , YM IN , XMAX , YMAX , ZN IN  ?  ZMAX  , DX  ?  0 Y  ?  DZ , NX » NY , NZ , N ; NORM ) 
200  STOP 
210  END 

220  SUBROUTINE  SCRSUP< SVECT, H,N) 

230  REAL  SVECT (3,7) >N(5) 

»  240  PRINT. ' INPIT  THE  NUMBER  OF  SOURCE  CURRENT  LOOPS  '5  INPUT. N 

250  1=1 

260  DO  10  WHILE  (I.LE.N) 

270  PRINT. 'INPUT  CURRENT  LOOP  SOURCE  '.I.'  POSITION  (XI,YI,ZI)» 

280  INPUT.  (SVECT(I.J) ,J=1»3) 

290  PRINT. 'INPUT  CURRENT  LOOP  WIDTH(X)  AND  LENGTH (Y) ' 

303  INPUT.  (SVECT(IiJ) . J=4.5) 

310  PRINT, 'INPUT  CURRENT  LOOP  '.I.'  STRENGTH (N*I ) » 

323  INFUT.il (I) 
iS  1=1+1 
34S  F1ZTURM 
joQ  END 

36»  SUBROUTINE  GRIDtXMIN.YMIN.XMAX.YMAX.ZMIN.ZMAX.NX.NY.NZ.DX.DY.DZ) 

270  REAL  XMIN.YMIN.ZMIN. XMAX ; YMAX . ZMAX , DX , DY , DZ 

38«i  INTEGER  NX.NY.NZ 

370  PRINT, 'FIELD  POINT  CRIO'JPRINT, 

403  PRINT, 'INPUT  MAX  FIELD  COORDINTES  (X.Y.Z) '5INPUT, XMAX, YMAX, ZMAX 
‘i13  PRINT, 'INPUT  MIH  FIELD  COORDINATES  (X.Y.Z)  ' 5INPUT, XMIN.YMIN.ZMIN 
423  PRINT. 'INPUT  GRID  RESOLUTION (DX.DY.DZ)  'HNPUT.DX.DY.DZ 
130  NX= (XMAX-XMIN) /DX+1 

: ‘0  ny=;ymax-ymin)/dy+i 

>'■=•3  NZ=(ZMAX-ZMIN)  'DZ+1 
460  END 

470  SUBROUTINE  SCRPRTi SVECT, M.N) 

480  REAL  SVECT (5 ,7) »M(5) 

%  490  INTEGER  N 

533  PRINT, 'SOURCE  CURRENT  LOOP  FIELD  INPUT  PARAMETERS' 

510  PRINT, .PRINT  100. 

520  130  FORMAT (2X, 'CURRENT  LOO?', 6X, 'XI', 8X»'YI',8X,'ZI',7X, 'WIDTH'. 

330  "5X, 'LENGTH' ,5X, 'N*I'/) 

S**0  PRINT  110,(I,{SVECT(I,vi)»J=l,5),M(I),I=l,N) 

-  350  110  FORMAT  CSX ,11, 5X » 5F 10. 4.F9 . 1 > 

560  RETURN 
573  END 

3 80  SUBROUTINE  GRBPRT (XMIN , YMIN, XMAX , YMAX , ZM IN, ZMAX , DX , DY , DZ) 

59?  REAL  XMIN, YMIN, XMAX, YMAX, ZMIN, ZMAX, OX, DY.DZ 
600  PRINT , 1 PRINT , 'FIELD  GRID  PARAMETERS' 

610  PRINT, 'XMIN  =  ' ? XMIN, 7  YMIN  =  MNIN.’ZMIN  =  ’iZMIN 
*  620  PRINT, 'XMAX  =  ' .XMAX, '  YMAX  =  ', YMAX, 'ZMAX  =  ',ZHAX 

63?  PRINT, 'DX  =  '.DX,'  DY  =  ',DY,'DZ  =  ',DZ 
640  RETURN 
65 0  END 

66 3  SUBROUTINE  DFCALC (SVECT ,M , XMIN, YMIN , XMAX , YMAX , ZMIN, ZMAX , DX , DY , DZ, NX , NY , 

570  «NZ,N,NORM) 

I  REAL  SVECT (5, 7), M(5), XMIN. YMIN, ZMIN, XMAX, YMAX, ZMAX, DX, BY, DZ, NORM 

.-.os  integer  nx.ny.nz.n 


703  REAL  X>Y>Z>SPX>$PYfSPZ>INTG>Iit  I2>I3>I4>X.JrYJT?.JtXL>YL»Bi4) 

710  CALL  SCRPRT(SVECTfHfN) 

723  CALL  GRBPRT i  XM IN  >  YMIN  r XMAX t YMAX » ZMIN » ZKftX ; DX >  BY » DZ) 

730  WRITE! 1 >500) 

743  2 =211  IN 

750  DO  US  L=liNZ 

763  Y=YHIN 

770  DO  100  1=1 > NY 

783  X=XHIN 

790  DO  90  J=1>NX 

833  S(ll=0?B(2)=0?B(3i=0 

810  BO  S0  K=1>N 

820  SPX=SVECT (Kr  1 ) ?XJ=X-SPX 

830  3PY=SVECT  (K»2)  >Y.J=Y-$PY 

840  SPZ=$VECT(K>3)>ZJ=Z-SPZ 

850  Xl=SVECT(K>4)/2 

363  YL=$VECT(K>5)/2 

870  I1=INTG(YJ>X.J>XL>ZJ>YL> 

833  I2=IH7G(XU>Y0rYLiZJ>XL) 

890  I3=INTC{YJ?XJ>”XL>ZJ>YL) 

900  I4=IN7G(XJiYJ>-YL>ZJ>XL) 

910  BU)=BU)+Z.J*(U-!3> 

923  B(2)=3<2)+ZJ*U2-I4> 

933  B(3)=3!3)-(XJ-XL)*Ii-{YJ-YL)*I2+(XJ+XL)'*I3+{YJ+YL)*I4 

940  30  CONTINUE 

950  BO  120  11=1 j3 

960  120  B ( I! ) =B ( I I ) /NORM 

970  B ( 4) =SQST (B(i) **2+B (2) **2+B (3)**2) 

980  IF (NORM. SO. 1.) 

990  NORM'S (4) 

1056  RETURN 
1010  END  IF 

1325  WRITE (1 >515)  XjYiZ.E 

1030  '500  FORMAT(r0*>5X>,Xf>10X>,YM0Xif2M3XtsBX,»9X»fBY'r9X>rBZf> 
1040  &9X>'B'/) 

1350  510  FORMAT t iX »3F5 0.2i 4Fi i . 3) 

1063  90  X=X+DX 
1970  WRITE U  >520) 

1080  100  Y=Y+0Y 
1090  WRITE'. i >520) 

1100  110  Z=Z+D2 
1110  520  FORMAT (/) 

1120  RETURN 
1120  END 

1140  REAL  FUNCTION  INTG(QrP>P0>ZiUL) 

1150  REAL  O>P>P0>2>UL>8»C>B 
1160  INTC=0 
1170  B=-2*Q 

1180  C=O«2+(P-P0)  **2+Z**2 
1196  D=B**2-4*C 
1200  IF(D.NE,2) 

1210  DO  10  I=2>3 

1 220  INTG=! NTG- (-1) **!*?* (2*UL+B ) /B/SGRT (UL «-2+B*UL+C ) 

1230  10  UL=-UL 
1240  ELSE 
1250  DO  20  I=2>3 

1263  INTG=INTG-(-l)**I*0.S5/(UL+B/2)**2 

1270  23  UL=-UL 

1280  END  IF 

1293  RETURN 

1300  END 
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COMPUTATION  OF  MAGNETIC  FIELD  LIMIT: 
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CRITERION;  V/V  >10 
A  — 

where  V  =  PLASMA  FLOW  VELOCITY 
VA  *  ALFVEN  SPEED 
“  B/(u0  MN)*5 

B  =*  MAGNETIC  FIELD  STRENGTH 

M  '  PLASMA  PARTICLE  MASS 

N  *  PLASMA  PARTICLE  DENSITY 

B  <  V  (y0  MNJ^/IO 

In  MRS  UNITS:  Wo  10“6 

M  -v  10-2 S 
N  1022 

V  -v  10* 

•  *•  B  1  io-2  (102  GAUSS) 
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